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F(z)=[F 0 (z 4 ) F,(z 4 ) P,(z 4 ) F 0 (z 4 )] 

3. Method according tojany of the claims 1 and^j characterized in that, 
in said sending filter (12), a filtering step followed by a step of interpolation by a 
factor of M=4 is performed. — — C^Ai^\ I - -~ 
5 4. Method according tojany of the claims 1 a 3~J characterized in that, in 

said reception filter (15), a step of decimation by a factor M=4 is performed, 
followed by a filtering step. — - Ol— ' ~~ "~ 

5. Method according to jany of the claims e 1 to 4^ characterized in that 
said sending filter (12) and/or said reception filter (15) have a structure in the form 

10 of at least one lattice. 

6. Method according to claim 5, characterized in that said sending filter 
(12) and said reception filter (15) are each constituted by a pair of polyphase 
components respectively given by the following equations : 

F J=gA(a JA(z)A(a„_ ,)... A(z)A(a 0 )[o] 

-F,l [01 
[ F J - gA(a n )A(z)A(a„_ , A(z)A(a 0 ) [,] 

f 1 a l ri o i 

* a) "L, and (2)= oz- 

with : L 1 J 

where g is a non-null constant of standardisation and otj are real coefficients. 

7. Method according to claim 6, characterized in that it implements a 
two-lattice structure. 

20 8. Method according to claim 6, characterized in that it implement a 

single-lattice structure working at a double frequency. 

9. Device for the filtering of Nyquist digital signals with null inter- 

s5'mbol interference designed to process a physical signal transmitted between a 

sender and a receiver through a transmission channel, 
25 based on an Nth order P(z) = F 2 (z) symmetrical filter implementing an 

oversampling factor M = 4 and forming a matched pair comprising a sending filter 

(12) and a reception filter (15), 



1 WBecmim UDECJOM 

lOp 09/719498 

Method for the making of digital Nyquist filters with null inter- 
symbol interference and corresponding filtering device 

The field of the invention is that of filtering of digital signals. More 
specifically, the invention relates to the making of Nyquist filters such as those 
used for example in baseband and single-carrier transmission. 

Nyquist filters play an essential role in transmission systems. In most 
cases, this filtering is distributed between sending and reception in the form of two 
filters known as Nyquist root filters. When sending, the aim especially is to limit 
the band of the signal sent. At reception, the filter must eliminate all the different 
noises that may damage the useful signal, especially the interference noise due to 
adjacent channels. In modern systems, these filters are made in digital form. 
There are many methods to synthesize these filters. A key element of this 
synthesis, apart from compliance with frequency specifications, is the need to 
obtain the smallest possible inter-symbol interference (ISI). 

To ensure the most efficient possible conditions of transmission, it is 
sought generally to obtain real coefficient filters which, ideally, will meet the 
following criteria: 

null inter-symbol interference; 
frequency selectivity; 

sending and reception filters forming matched pairs; 
phase linearity. 

Appendix 1 describes these different aspects and details and the corresponding 

constraints. This appendix, just like the following ones, is of course an integral 

part of the present description. 

A complete transmission system must generally also carry out other 

indispensable functions such as: 

analog filtering at sending which is designed to reject the output 
harmonics from the digital filter as well as analog filtering at 
reception whose function is to limit the spectrum of the signal to be 
sampled solely to the useful band; 

digital-analog conversions (DAC) at sending, and analog-digital 
conversions (ADC) at reception. 
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The invention essentially relates to the obtaining of digital filters, and no 
detailed description shall be given of the functions relating to the analog part. It 
must be noted however that, in practice, the analog filtering to be implemented is 
substantially less restrictive in terms of steepness of the filters than digital filtering. 
5 More specifically, the complexity of the embodiment of the rejection filter at 

sending is directly related to the value of the oversampling factor chosen. 

As for the converters, the most notable fact from the viewpoint of filtering 
is that of digital-analog conversion which introduces a sine(x)/x filtering that may 
be conventionally compensated for by the analog filtering. This processing can 
10 also be done digitally. 

Furthermore, for reasons of implementation, the filters are generally 
sampled at double or quadruple frequency. 

Figure 1 illustrates the digital part of a communications system of this 
kind. The signal to be sent X(z) 11 supplies the sender 12 which delivers a 
1 5 filtered signal Y(z) 13 transmitted through a transmission channel 14 to a receiver 

15 that delivers the output signal S(z) 16. The signal X(z) 11 is first of all 
subjected to oversampling 121 and then to a filtering 122 F t (z). In the receiver, 
the signal undergoes a reception filtering F r (z) 151 and then a decimation 152. 

With T as the elementary delay linked to the sending and reception digital filters, 
20 we have T s = MT, where T s is the symbol duration and M the factor of 

oversampling. 

The value of having a high oversampling factor is especially the fact that it 
facilitates the making of the analog filter that will follow. One drawback is that it 
correspondingly increases the speed of the converter (DAC). In practice, the most 
25 reasonable compromise corresponds to a value of M = 2 or 4. The invention 

relates specifically to the case where the oversampling factor is M = 4. 

Nyquist filter construction techniques are already known. However, the 
different classes of known filters always entail the release of one of the following 
desired properties : zero inter-symbol interference, phase linearity (at least at 
3 0 sending) and the pair of matched filters. 

Often, in communication systems, it becomes necessary to accept that the 
filtering introduces a non-zero ISI. This however sometimes raises problems, for 
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example in the case of modulations with very large numbers of states (64, 256, 
...). 

It is an object of the invention especially to overcome these drawbacks of 
the prior art. 

5 More specifically, it is an object of the invention to provide a method for 

the making of Nyquist filters with null inter-symbol interference that also meets 
the desired conditions of frequency selectivity and phase linearity in a matched pair 
configuration. 

Another aim of the invention is to provide a method such as this that 
1 0 enables easy practical implementation of the filters obtained and takes account of 

certain practical constraints such as the effects of the quantification of the 
coefficients. More specifically, a goal of the invention is to provide a method to 
ensure the maintenance of the property of cancellation of inter-symbol interference 
after this quantification. 

1 5 Another aim of the invention is to provide a method of this kind that can be 

implemented by means of a cascade structure, especially in lattice form. 

These aims as well as others that shall appear hereinafter are achieved 
according to the invention by means of a method for the making of a digital 
Nyquist filter with null inter-symbol interference designed to process a physical 

2 0 signal transmitted between a sender and a receiver through a transmission channel. 

This filter is an Nth order P(z) = F^(z) symmetrical filter implementing an 
oversampling factor M=4 and forming a matched pair comprising a sending filter 
and a reception filter whose polyphase breakdown of F(z) can be written as 
follows: 

25 F(z) = F 0 (z 4 ) + z- 1 F 1 (z 4 )+ z- 2 F 2 (z 4 )+ z" 3 F 3 (z 4 ). 

According to the invention, N is chosen to be different from 4n, n being an 
integer and the coefficients of the polyphase breakdown of F(z) are such that: 
IfN=4n+l, F 1 (z)F 1 (z)+ z" 1 F,(z)F,(z) = yz~" 
If N=4n+2, 2F 0 (z)F 0 (z) + F°" (z) + z~ 1 F 2 3 (z) = yz - " 
30 . IfN=4n+3, F 0 (z)F 0 (z) + F I (z)F ! (z) = yz" n 

F being the mirror symmetry of F and y being a non-null constant. 
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Thus, the method of the invention, by construction, ensures that the ISI is 
perfectly null. 

Preferably, N is equal to 4n+3 or 4n+l and: 

said sending filter performs an interpolation by a factor M = 4 and has a 
circuit arrangement corresponding to a polyphase breakdown known as the 
type II breakdown, such that: 



F(z)=[z" 3 z 



Fofe ) 
Fife 4 ) 
F.fz 4 ) 
F 0 (z\ 
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and said reception filter performs a decimation by a factor M = 4 and has a 
circuit arrangement corresponding to a polyphase breakdown known as the 
type I breakdown, such that: 



15 



20 



F(z)=[F 0 (z 4 )F,(z 4 ) f,(z 4 )F 0 (z 4 )] 



The filters thus obtained are simple to make and implement. 

Advantageously, in said sending filter, there is performed a filtering step 
followed by a step of interpolation by a factor M = 4. Similarly, in said reception 
filter, advantageously a step of decimation by a factor M = 4 is performed, 
followed by a filtering step. 

This construction (permutated structure) reduces the rate of the operations 
by a factor 4. 

According to a preferred embodiment of the invention, said sending filter 
and/or said reception filter has a structure in the form of at least one lattice. 

Indeed, in this form of circuit arrangement, the constraint of perfect 
reconstruction is structurally integrated. 



Advantageously, said sending filter and said reception filter are each 
constituted by a pair of polyphase components respectively given by the following 
equations: 



where g is a non-null constant of standardization and are real coefficients. 

According to one particular embodiment of the invention, the method 
implements a two-lattice structure. According to another embodiment, it may 
implement a single-lattice structure working at a double frequency. 

The invention also relates of course to filtering devices obtained by means 
of the above-described method. These are therefore devices for the filtering of 
Nyquist digital signals with null inter-symbol interference designed to process a 
physical signal transmitted between a sender and a receiver through a transmission 
channel, based on an Nth order P(z) = F 2 (z) symmetrical filter implementing an 
oversampling factor M = 4 and forming a matched pair comprising a sending filter 
and a reception filter, the polyphase breakdown of F(z) of this symmetrical filter 
being written as follows: 



F(z) = F 0 (z 4 ) + z- 1 F 1 (z 4 )+ z~ 2 F 2 (z 4 )+ z" 3 F 3 (z 4 ). 
According to the invention, N is different from 4n, n being an integer, and: 



Other features and advantages of the invention shall appear more clearly from the 
following description of a preferred embodiment given by way of a simple, 
illustrative and non-restrictive example and from the appended drawings, of 
which: 



with : 





Figure 1 , already commented upon in the introduction, provides a 
diagrammatic illustration of the digital part of a communications 
system; 

Figure 2, commented on in Appendix 1, represents a typical 
frequency specification of a sending filter; 

Figure 3, commented upon in Appendix 2, gives a general view of 
a bank of filters with two sub-bands; 

Figure 4 illustrates an embodiment of the analysis part of a bank of 
para-unitary filters with two sub-bands of Figure 3 in lattice form: 
Figure 5 recalls the structure of a sending/reception filtering system 
in the linear and para-unity case commented on hereinafter; 
Figures 6a and 6b give a view, for a 4n+3 order filter of two forms 
of circuit arrangement of the sending filter of Figure 5, using a 
polyphase breakdown respectively in a direct structure and a 
permutated structure; 

Figures 7a and 7b provide an illustration, in the case of the 4n+3 
order, of the reception filter of Figure 5 also in the context of 
polyphase breakdown respectively according to a direct structure 
and a permutated structure; 

Figures 8a and 8b show a lattice block of the sending/reception 
filters of Figures 6a, 6b, 7a and 7b respectively in the form of a 
direct lattice and a lattice with inverted outputs; 
Figure 9 gives a detailed view of a 4n + 3 order Nyquist root 
sending filter with two lattices implementing the inverted lattices of 
Figure 8b; 

Figure 10 shows a 4n+3 order Nyquist reception filter 
implementing four lattices; 

Figure 10b shows a reception filter implementing a direct lattice 
according to Figure 8a (see Appendix 4); 

Figure 11 recalls the equations used in the case of the inverted 
lattice of Figure 9; 



Figure 12 shows another embodiment of the sending filter 

implementing a single lattice and working at a double rate; 

Figure 13 shows an exemplary response in optimized frequencies 

of a 43rd order Nyquist root filter according to the invention; 

Figure 14 illustrates the embodiment of a block Z in the case of a 

4n+2 order filter as discussed in Appendix 3; 

Figure 15 shows the circuit arrangement block of the matrix 

also in the context of the embodiment of Appendix 3; 

Figure 16 shows the full circuit arrangement diagram of a 4-input 

and 4-output block for the 4n+2 order referenced Ma and discussed 

in Appendix 3; 

Figure 17 shows an embodiment of a half Nyquist sending filter 
for the 4n+2 order; 

Figures 18 and 19 (block N of Figure 18) illustrate a half Nyquist 
reception embodiment for the 4n+2 order discussed in Appendix 3; 
Figures 20 and 21 respectively show the sending and reception 
Nyquist root filters for the 4n+l order (or 4n+5 order) as 
discussed in Appendix 5; 

Figures 22 and 23 illustrate the results obtained by means of the 
mode of generic synthesis commented upon in Appendix 6. 
As mentioned already, the method of the invention can be used to make 
filters with M = 4 verifying the totality of the following criteria: 
ISI null by construction; 
linear phase filters; 

sending and reception filters forming a matched pair. 
Appendix 2, after developing the analogy between a Nyquist pair with M = 
4 and a bank of filters orthogonal to two sub-bands; shows that the two situations 
may occur according to the order of the filter. More specifically, it is shown that 
there is no solution for the N = 4n order filters and that the cases N = 4n + 1 and 
N = 4n + 3 are favored (theorem 2) in that they produce similar polyphase 
components. The condition for obtaining null inter-symbol interference is 
specified in the theorem 1. 
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The case of the 4n + 3 filters is described hereinafter. It proves to be the 
case used to obtain the simplest embodiment. The case 4n + 2 is processed in 
Appendix 3 and the case 4n + 1 (or 4n + 5) in Appendix 5. 

As a rule, the basic diagram implemented in shown in Figure 5. As 
5 compared with the initial general diagram of Figure 1 , it may be noted that Fp(z) 

122 and Fj^(z) 151 are replaced by F(z). This is indeed a linear phase matched 
pair. 

We must also note the introduction of a delay line z-r51 where r is a 
positive integer value dependent on the order, which in all cases makes it possible 
10 to obtain a null ISI with a given delay. The multiplier factor g 52 is used to ensure 

that, except for the delay, the output is quite identical to the input. 

In the case where N = 4n + 3, because of the relationships between the 
different polyphase components F(z), the expression (20) of Appendix 2 can be 
rewritten as follows: 
15; F(z) = F 0 (z 4 ) + z~ l F l (z A ) + z' 2 ^) + r 3 F 0 (z 4 ) (31) 

It is known that, apart from the delay, the sending and reception filters are 
identical. However, if it is sought to draw the best advantage of the multi-rate 
processing operations, namely interpolation at sending and decimation at 
; reception, the structural drawings will not be exactly the same. 

20 At sending, where an interpolation by a factor M = 4 is made, the forms of 

circuit arrangement shown in Figures 6a and 6b correspond to those known as 
type II polyphase breakdown [4] (the references cited are assembled in Appendix 
7), given in the case of the invention by: 



F(z)=lz" z" z" 1. 



(32) 



F L (z ) 

25 i F °^). 

In one of the depictions (Figure 6b), it is possible to note the permutation 
of the expansion (61) and filtering (62) operations used to carry out an equivalent 
processing by working at a rate four times lower than in the case of the initial 
diagram (Figure 6a). 
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In the reception part, the filtering is associated with a decimation operation 
for which a more appropriate writing of F(z) is that of what is called the type I 
[4] polyphase decomposition which in the present case can be written as follows: 



F(z)-[F 0 (z 4 )F,(z 4 ) F^FJz 4 )] 



(33) 



It will be noted that, with reference to Figure 5, the time limit is such that r 
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15 



Two variants of the circuit arrangement of this equation are given in 
Figures 7a and 7b (symmetrically with Figures 6a and 6b). 

The analogy between the 4n+3 order filter F(z) and the 2n +1 order bank 
of filters described by the theorem 3 (Appendix 2) is used also to rewrite the 
equation in the form: 

F n (z) -F ,(z)" 

= gA(a ^ ) A(z)A(a ^ ^...A (z)A(a q ) 

(34) 

The two pairs of polyphase components which are mirrors of each other are 
therefore deduced from each other by the following two equations: 



F(z) F (z) 



L° = gA(a „)A(z) A(a„_ , ). . .A(z)A(a 0 ) [J 1 



(35) 



. j -gA(a n )A(z)A(a n . ] )...A(z)A(a 0 ) 

l F oJ m (36) 

In each of these cases, the lattice structure shown in Figures 8a and 8b is the 



The two variants, with direct or inverted outputs of Figure 8, are then used to 
set up the sending and reception filters according to the diagrams respectively 
shown in Figures 9 and 10. It will be noted that the symbols at output of the 
filters indicate only which are the active filters. 
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By using the polyphase decomposition of the input signal, it is then possible 
to determine an equivalent diagram with only one lattice which however, in this 
case, is constrained to work at the double frequency. 

In the diagram proposed, the inverted lattice is presented, it being known of 
course that the direct lattice leads to a diagram of the same type. Figure 1 1 recalls 
the basic equations of this lattice. 

The output signal of the sending filter is written as follows: 

Y(z) = z _1 [F 0 (z 4 ) + z~ 1 F 1 (z 4 ) + z- 2 ^a 4 ) + z" 3 F 0 (z 4 )]x(z 4 ) (37)> 
which can be rewritten as follows: 

Y(z) = z' ] [F () (z 4 ) + z- 2 F i (z 4 )}X(z")+ Z - 2 \F x {z 4 ) + z- 2 F 0 (z*)]x(z 4 ) 

= [F 0 (r) + z-^( 2 2 )]z-' 

A division by 2 of the expansion factor (T 4 = T 2 T 2), and then a 
permutation of the expansion of 2 with the filtering produces filtered outputs such 
that: 

Y 0 (z) = [F 0 (z 2 )X (z 2 ) + z- 1 F , (z 2 )X (z 2 ) ] 
Y , (z) = [F , (z 2 )X (z 2 ) + z - 1 F 0 (z 2 )X (z 2 ) ] 

(39). 

In Yi(z), it is possible to recognize the higher output of the inverted lattice 

working at double rate and, plus or minus the sign of the term as a factor of z"*, in 
Y Q (z) it is possible to recognize the lower output. To recover exactly the same 

expression, it is therefore enough to multiply this output of the inverted lattice by 
(-D n 

One embodiment based on this principle is shown in Figure 12. 
For the reception part, it is possible to consider two cases: 
- the general case, with four lattices, which will need all the desired 
properties (Figure 10); 

- the simplified diagram for which only the null ISI and the phase linearity at 
sending are guaranteed (Figure 10b). 

This particular aspect is discussed in Appendix 4. 

In the case of Figure 9, the system works at the lowest frequency. The cost 
of each filter corresponds to a lattice structure with n+1 cells plus that of a sign 
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inversion. To obtain a unit gain, a multiplication by the factor g generally must be 
accounted for in the diagram of the reception filter. 

The computations for each elementary filter correspond to the two multipliers 
and adders. However, as in the case of each lattice, one of the inputs is at zero, 
5 and a multiplier and an adder can therefore be eliminated. The cost is therefore 

one multiplier and one adder for this first cell and double for the following, 
namely (2n + 1) multipliers and (2n +1) adders. For each filter, we will therefore 
have (4n + 2) multipliers and (4n + 2) adders. At the lowest frequency, we will 
therefore have to make (4n + 2) MPU and (4n +2) APU. 
1 0 For the reception part, in the case of the simplified diagram illustrated by 

Figure 10b, the operational complexity is that of a single lattice rather than that of 
three adders, the multiplier by g and the inverter. The computations performed at 
the lowest rate can be estimated at (2n + 3) MPU, (2n + 5) APU and one 
inversion. 

15 In the general case illustrated by Figure 10, at the lowest frequency, the 

complexity is (8n + 4) MPU and APU and one inversion. 

The method of synthesis is performed for example in three steps, illustrated 
here below by the computation of a 43rd order filter. 

Step 1: Using the synthesis tools [3], it is possible to obtain a first 

20 set of transversal coefficients that are close to null IS I. 

These are the initial coefficients of Table 1. 
Step 2: B 4 y identification, we can compute the corresponding initial 

lattice coefficients given in the Table 2. Since the initial set 
of transversal coefficients is not exactly at null ISI, a new 
25 set of transversal coefficients is thus obtained known as 

reconstructed coefficients (cf. Table 1). 
Step 3: A local optimizing of the initial lattice coefficients then 

provides for slight improvement in the frequency response 
of the filter. Tables 1 and 2 give the values of these lattice 
3 0 and transversal coefficients after this optimization (Tables 1 

and 2). 



12 



Table 1 : Transversal coefficients of the 43 rd degree symmetrical filters 



10 



15 



% 


Initial 


Reconstructed 


Optirnizal 


0 


0.006302365568 


0.006302365568 


0.00374315552336009 


1 


-0.01335763186 


-0.01335763186 


-0.00764365688411938 


2 


-0.003586256877 


0.000139948S24594686 


0.000256865997947587 


3 


0.003281590529 


0.000296616382849522 


0.00052452951560661 


4 


0.008254154585 


0.00825445268157875 


0.00210146572675043 


5 


0.008753151633 


0.00875251982827294 


0.00961384935212083 


6 


0.002954708412 


-0.000818038274240153 


-0.000565430770503378 


7 


-0.006948650815 


-0.00231664836261536 


-0.00210883853472157 


8 


-0.01474214438 


-0.0147620230725521 


-0.00994253905900314 


9 


-0.01369718742 


-0.0136538137768085 


'-0.0139493173848354 


10 


-0.001682000235 


0.00199238622180159 


0.00140539009490696 


11 


0.01571460254 


0.00939101589973353 


0.00785652590927169 


12 


0.026870884 


0.0269017853424573 


0.0234276140713132 


13 


0.02112355083 


0.0209715158594175 


0.0216783222304717 


14 


-0.003340088297 


-0.0032097877874163 


-0.00339966101869637 


15 


-0.03511243314 


-0.0308000965511264 


-0.0274417538899884 


16 


-0.05300379917 


-0.0530202725288858 


-0.0511658629696837 


17 


-0.0364661812 


-0.0360408764445379 


-0.0377795925012697 


IS 


0.022169S768 


0.0140041435375816 


0.0161657062857809 


19 


0.1105730906 


0.116698398397514 


0.109490464868349 


20 


0.1994492114 


0.199658215352476 


0.199471282748299 


21 


0.2550483644 


0.253344105750997 


0.257329314844394 



Table 2 Coefficients of the lattice structures 





Initial 


Optimizal 




0.006302365568 


0.003743155523 


cti 


2.1194631945539 


2.04203561311232 




- 0.7583059650965 


- 0.71854479175763 


&3 


0.6811397085317 


0.73802350460S93 


a 4 


- 0.4637577935126 


- 0.55360950120952 


&5 


0.4470434662649 


0.51415010432349 


a< 6 


- 3.7431288142180 


- 5.34790460147966 




- 1.5196264288798 


- 1.52189568702330 


ct 8 


1.2422617481334 


1.44910973460487 


a- 9 


- 0.5612421724760 


- 0.63596605846521 




0.15S8040486787 


0.18869460160260 


Qlj 


- 0.0222057611677 


- 0.068622S4945003 
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The frequency response for the 43rd-order filter computed in this way is 
shown in Figure 13. 

This is a first result that may be further refined. Furthermore, of course, 
other computation techniques may be considered for different technical orders and 
constraints. 

It will be recalled especially that the case N = 4n + 2 is presented in the 
Appendix 3. Furthermore, Appendix 5 discusses the case N = 4n + 1 (or 4n + 5). 

Details of the assessment of the complexity of the reception filters and a 
method of generic synthesis are developed in Appendix 6. It will be noted that 
this method of synthesis can be applied to all three types of approaches described 
here above. 
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APPENDIX 1 

The distribution of the Nyquist filtering on two filters therefore commonly called 
Nyquist root filters or, in short, half-Nyquist filters, must take account of the constraints 
related to the ISI, the frequency specifications, the linearity of the phase and what is called the 
5 matched pair property. In this appendix, we define these different notions assuming, to 
simplify the notation, that T- 1. Furthermore, in the context of the embodiment, we shall deal 
only with the case of the finite pulse response filters (FIR). 

1 - Intersymbol interference (ISI) 

10 Let x(n) be the input signal and X(z) its z transform, namely X( z ) =X„ M n )z' n and let 

it be agreed that all the signals of the transmission system will be referenced according to the 
same principle. The null ISI constraint amounts to the dictating, when there is no noise, of the 
output signal referenced S(z) which is identical, except for a delay, to X(z). 

Let Y(z) be the output signal from the sender. It can be written as follows: 

15 

= Y{z) = F 7 (z)X{z M ), (1) 

At output of the reception filter, the signal referenced U(z) has the following expression: 
20 U{z) = F T {z)F R (z)X{z M ). (2) 

The output signal is therefore quite simply: 

i M ~ i , - r i M ~ i , , i 

= Yl f Q U ^ w ^ = [Jf f Q Fr(^w k )F R (z^w k ^ X(z), (3) 
2 jn 

where w = e is the Mth root of unity. The null ISI condition can then be summarized in 
the following equation: 

1i;f t ( : »,')f ! ( ! %') = .--', (4) 

where d corresponds to the delay introduced by the two filters. 

Let P(z) be the filter produced, i.e. P(z) = Fj(z)F R (z). In a formal method commonly 
used in with multi-rate signal processing, the equation (4) can then be written as follows: 



25 



30 
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P{z)i M =z~ d . (5) 

A filter P(z) that verifies this relationship is a Nyquist filter or again an Mth band filter. 
5 Let n T and n R be the respective orders of the sending and reception FIR filters. The 

filter P is therefore an n P = n T + n R order filter and can be put in the following form: 

np 

^) = EK»)~~- a - (6) 

10 If the equation (5) is not verified, it is commonly the "distance" that is verified with 

respect to the property by either of the following expressions: 



1 5 where 



(7) 



f n 2 = -±- ZAd-kM). (8) 

20 2. The frequency specifications 

Figure 2 shows a typical frequency specification for digital transmission filters: 

• The roll-off factor defines, for these lowpass filters, the passband by 

TL TZ 

[0, a> p = jj- (1 - p)} and the attenuated band by [co s = j^- ( 1 + p), n] 

• The passband ripple specifications referenced 8 } and an attenuated band referenced 8 2 
25 are generally the same for the sending and reception filters. 

• At the pulsation n/M, generally the constraint F R (n/M) = F^n/M) = is 
dictated. 



3 - The phase linearity 

3 0 The choice of characteristic phase frequency filters that are perfectly linear is generally 

recommended for digital transmission systems [1, p. 325]. This especially has the advantage 
of preserving the instants of passage through zero of the transmitted binary trains. 
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4 - The matched pair of filters 

In the case of linear modulation, if we assume that ISI is zero for the entire transmission 
system and that the channel noise is additive, white and Gaussian (BBAG), it is known that the 
matched pair, namely Fj(z) = z~ N F R (z~ l ) with N = n T = n R , the order of each filter, is optimal 
for the criterion of the signal-to-noise ratio (SNR) [2, pp. 51 sq.]. The SNR which is equal to 
~T\ \f iPl'F IP t * ien reacnes i ts maximum with ||^Rli 2 ||-Pr|p = 1. 
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APPENDIX 2 

1 - Orthogonal filter banks with two sub-bands 

A filter bank with two sub-bands can be shown according to the diagram of Figure 3. 
In order to meet the condition of non-aliasing, the synthesis filters G 0 (z) and Gj(z) are directly 
deduced from the analysis filters by G 0 (z) = Hj(-z) and G } (z) = -H 0 (-z). 

It can then be shown [4] that the condition of perfect reconstruction (PR) can be 
expressed solely from the polyphase matrix H[h 0 ,hJ (z) of the analysis bank: 

where the values of H kl (z) are the polyphase components for each filter H k (z), namely: 
H t (z) = £z->H kJ (z 2 ) m 

1=0 

In the case of the orthogonal banks with two sub-bands, the filters are 2n+l odd-order 
filters and the polyphase components are all of the same n degree. The bank is then PR 
if the determinant of the polyphase matrix is a monome. In other words if: 

20 Det H [FoJ7l] (~) = H oja (z)H u {z) - H lfi (z)H 0 .i(z) = ^ (11) 

with P as a non-null constant. 

Let Q (z) be the mirror filter of the mth order filter Q(s), namely: 
Q (s) = z- m Q(z- ! ). The orthogonal bank is characterized by filters H 0 (z) and H } (z) said to be 
25 quadrature conjugates, which can be written as follows H,(z) = ±z- (2n+1) H 0 (z' 1 '). Hereinafter, 
unlike the presentation given in [4], we choose the writing with a plus sign. Thus, the 
polyphase components meet the following relationship: 



15 
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Fcu(^) = -#i,o(~) , #i.i(~) = #o.o(~). (12) 
Consequently, the polyphase matrix can be written as follows: 
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and the relationship (11) can be written in the form: 

Det H [Ho , ffl] (c) = H 0fi {z)H 0fi (z) + S lfi (z)H lfi (z) = /fcr" , (14) 

The output signal is then such that X ( s ) - yz- (2n+1) X(z) with y as a non-null constant. 
The matrix H[h 0 ,h,] (z) is para-unitary and can be cascaded with [4]: 

H [ffoJ?l3 (-) - 5 A(Q n )A(^)A(a n _ I ) . . . A(z)A(a 0 ). (15) 

where g is a non-null standardization constant, Oq, a n are n+1 real numbers and where the 
matrices A(oc) for a real and A(z) are defined by: 

M a )=[l a °] ■ A ( .-)=[j (16) 

15 

An embodiment of a synthesis bank based on this mathematical writing produces the 
lattice depiction of Figure 4. 

2. General characteristics of the case of the linear phase matched pair for 
20 M=4 

For the linear phase matched pair, we necessarily have Fj(z) and F R (z) symmetrical 
such that Fj(z) = F R (z), which can be written with only one expression F(z). Let Abe the 
order of F(z). We therefore have: 

25 F0O = I>s-* (17) 

fc=0 

The produced filter P(z) is then expressed by: 

pw = ^) = E^ (18) 

with of course p(k) = p(2N - k). P(z) meets the zero ISI condition (5) if: 



p(k) = 0foik-N=4l,l±0 



(19) 



19 
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It is possible to make a first observation with regard to the N - 4n order filters. 
Proposition 1 - Let F(z) be a 4n order symmetrical filter. F(z) cannot be a null 
ISI filter. 

In order that F(z), a symmetrical filter, may be a null ISI filter, we must have/ 0 ?0. 
The highest degree monome of P(z) is fo~ z~ 8n while its 4n degree central term is z A . Since the 
difference in degrees is a multiple of 4, F(z) is not at null ISI. 

The polyphase decomposition of F(z) for M = 4 can be written in the following form: 



F(=) = F 0 (z 4 ) + z-'-F^z 4 ) + z~ 2 F 2 (z 4 ) + z^Fsiz 4 ) . 



(20) 



(21) 



To carry out the analysis of the different possible cases as a function of the Nth order, 
the property of symmetry of F(z) and the fact that 4nj + l<N will be used. 

• If N = 4n, the degree of F 0 (z) is equal to n and that of the values F^z), i = 1,2,3 is 
strictly below n. According to Proposition 1, it is known that this case is of no use 

2 0 whatsoever for our problem. 

• If N = 4n + 1, Fj(z) is an n degree value, F 0 (z) is a < n degree value while F 2 (z) and 
F 3 (z) are < n - 1 degree values. Fj(z) then corresponds by mirror symmetry to F 0 (z). Fj(z) 
is therefore actually an n degree value. We therefore have: 

25 F 0 (z) = et F 3 (z) = F 2 (~~) (22) 

• If N = 4n + 2, we get n 2 = n, n x = n and n 3 = n- 1. F 2 (z) is symmetrical to F 0 (z). 
We therefore have: 

30 F 2 (z) = F 0 {z) , F^^F^z) , F 3 (=) = F 9 {=) (23) 



♦ If N = 4n + 3, 4n, + 1 < 4n + 3, VI (0 < / < 3), and the filters F t are all < n degree 
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filters. Indeed, the highest-degree term of F 0 (z) is the 4n + 3 degree term giving an n degree 
term at F 3 (z). The highest-degree term of F 3 (z) corresponds by symmetry to the constant term 
of F 0 (z). We therefore have: 

5 F,(z) = F Q (z) et F 2 (z) = F l (z) (24) 

Note: In the equation (22), we have deliberately written F^z) = F } ( z ) and not F } (z) = 
F 0 ( z) because it is certain that the degree of Fj(z) is n since the constant term of F^z) is not 
zero but not that the degree of F 0 (z) is n. Furthermore, this case where f 1 =f N .j = 0 is 
1 0 considered further below. 

The product F 2 (z) is then developed in taking account of the decomposition (20) and the 
following characterization of the filters F(z) verifying the null ISI property is obtained. 

Theorem 1 - Let F(z) be an Nth order symmetrical filter verifying F(z = 0) 0. 
Then, depending on the values of N, F(z) is at null ISI if and only if: 
15 • N = 4n + 1, 

F 1 (z)F 1 {z) + z-'F 2 (z)F 2 (z) = 72 r» , (25) 
• If N=4n + 2, 

20 

2F 0 (z)Fo(z) + F;(_~) + z-'Fi(z) = , (26) 



• If AT=4n + 3, 
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Fo(z)F 0 (z) + F 1 (z)F 1 (z) = 1 z~ 



(27) 



where y is a non-null constant. 

The following theorem shows that all the N = 4(n - 1 ) + 3 order solutions are obtained 
from N = 4n + 1 order solutions. 
3 0 Theorem 2 - Let F(z) be a null ISI symmetrical filter of the 4n + 1, n >1 order. 

Then ifF/z), i = 0, 3 designates the polyphase components, F 0 (z) is of a degree strictly 
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below n and F } {z) can be written in the following form: Fj(z) = z^K^z). We then have 

K } ( z ) zz F 0 (z) and the filter F (z) with polyphase components [F 0 (z), F 2 (z), F 3 (z), K } (z)] is 
a symmetrical 4(n - 1 ) + 3 order, null ISI filter. 

Demonstration - This immediately results from the relationship (25) since the 
5 constant term of the member on the left is null. We therefore have 

F/z = 0) F j( z - 0) ^fofj = 0. Since Fj(z) contributes with/ 0 to the highest-degree 
term,/ 0 ^0 and therefore/; = 0. Writing Fj(z) in the form Fj(z) = z-'Krfz) and 
knowing that F j(z) = K } (z), the relationship (26) can also be written as follows: 

1 0 z-'K^K.iz) + z- 1 F 2 (z)F 2 (z) = T _~- , (28) 

Since & } (z) = F 0 (z) and therefore K/z) = ^ 0 (z), we obtain the relationship (27) for a 
filter where Kj(z) is substituted for F 3 (z) and therefore has the polyphase components: [F 0 (z), 
F 2 (z), F 3 (z), Kj(z)]. It is therefore possible, in one example, to verify the passage from a 9th 

1 5. order filter to a 7 th order filter. 

Example 

Fo(z) =f 0 + f4Z~ 1 +f 1 Z- 2 ^fo + f.Z- 1 

F^z) = + fa- 1 + fa' 2 =2- 1 (/4 + /0^ 1 ) =Z~ 1 K 1 {z) 

Q F 2 (z) =/ 2 + / 3 ^ ( 29 ) 

F 3 (z) ^fs + fiz- 1 

2 0 The composite filter of [F 0 (z), F 2 (z), F 3 (z), Kj(z)] is a 7th order filter and is then 

expressed by: 

fo + I2Z- 1 + /a-"" 2 + + hz'"- + hz~'> + .hz' e + foZ~ 7 (3Q) 



25 A third theorem shows that the N = 4n + 3 order solutions are obtained from the 

orthogonal two sub-band filter bands. 

Theorem 3 - Let [H 0 (z), H } (z)J be a bank of2n + 1 order orthogonal filters with 

two sub-bands whose polyphase matrix is given by (13). Then the filter H (z) for which 
the four polyphase components are [H 0 0 (z), H 10 (z), H l0 (z), H 0 ,o(z)] is a linear phase, 4n 
30 + 3 order, null ISI filter. Reciprocally, all the filters of this type are obtained from 
orthogonal filter banks with two sub-bands. 
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Demonstration - The demonstration is immediate since the 4n + 3 order linear phase, 
null ISI filters are characterized by the relationship (27) which, in the terms of the theorem, 
becomes equivalent to (14). 
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APPENDIX 3 
N - 4n + 2 order Nyquist pair 

1 - Initial factorization 

Theorem 1 in Appendix 2 gives us, inter alia , the condition that a 4n + 2 order F(z) 
filter must meet in order to be at null ISI. It may be recalled that it is written as follows: 



2F 0 (z)F 0 (z) + F*(z) + z -i F fcy_ 



(31) 



where y is a non-null constant and where the components F t (z) are all nth degree components 
apart from F 3 ( z ), and meet the relationships F 2 (z) = F 0 (z), fJz) = F x (z) and F z (z)= F 3 (z). 

Theorem 4 - Let a t , i = 1, n + 1, n > 0, n + 1 be constant values. We consider 
the column vector [F 0 ( z ), F/z), F 2 (z), F 3 ( Z )] T defined by the following equality: 

Friz) . 

F 2 (z 



Fziz) 



-- M{a x ) Z M(a 2 ) Z ... Z M(a n ) Z 



(32) 



Then the filter F(z) whose polyphase components are the values Ffz), i = 1, 4, is a null 
ISI, 4n +2 order, monic, linear phase filter. 

Reciprocally, any null ISI, linear phase, monic, 4n + 2 order filter accepts a 
decomposition of this form. 

Demonstration - For n = 0, the polyphase components are the constant polynomials 
F 0 (z) = 1, Fj( z ) = a 1? F 2 (z) - 1 and F 3 (z) = 0. These polynomials trivially verify the 
condition (31) as well as the conditions of symmetry and degree. We then have F(z) = 1 + 
ajZ A + z 1 . Let us now suppose that the components of a 4n + 2 order F(z) filter F t (z), i = 0, 
3 verify the condition (31) as well as the conditions of symmetry and degree. Let us build 
L ; (z), i = 0, 3 with: 



Lo(z) 

L 2 (z) 

Ls(z) 



= M( ai ) Z 



Fo(z) 
F.{z) 
F*(z) 
Fs(z) 



(33) 



L 0 (z) = F 0 (z) - 



''F^ + a^Fsiz) 



(34) 



24 



Li(z) = ~ ai F 0 (z) - a^F 2 {z) + (1 - ^)z~ l F 3 {z) , (35) 

L 2 (z) = -^Fo(z) + z~ 1 F 2 {z) + a.z^Fsiz) , " (36) 

Ls(z) = (1 + ^)Fi(s) - ' (3?) 

The verification of the conditions of degree of symmetry on the components L t (z), i = 
0, 3 is done almost immediately. By developing the analogous expression of (31) for the 
components Liz), i = 0, 3, we obtain: 

10 2L 0 (z)Lo(z) + Ll(z) + z^Ll(z) = 

(2Fo(z)F 0 (z) + Fftz) + z-'Fi(z)) (38) 

According to (31), the direct part of the theorem is therefore demonstrated. 

The reciprocal is demonstrated if it can be established that the relationships (34-37) can 

1 5 uniquely determine a coefficient a } and the components Flz), i - 0, 3 of a null ISI, linear 
phase, monic F(z) filter. Since we should have F 0 (z = 0) = 1 (since F(z) is monic and 
symmetrical, its constant term is equal to 1), according to (35) we have a 2 - -L : {z - 0). If a ; 
= 0, then we are in the particular case where L 0 (z) is an n degree term, the n + 1 degree term 
being zero, the symmetiical polynomial L t (z) has a null constant term and its highest degree 

20 term is null, and finally the constant term of L 2 (z) is null (but not its n + 1 degree term whose 
coefficient is 1). The relationships (34-37) are then inverted by F 0 (z) = L 0 (z), F } (z) = L 3 (z), 
F 2 (z) = L 2 (z) quo z- 1 and F 3 (z) - L } (z) quo z' 1 . It is now assumed that a 1 # 0 but this is not 
actually necessary. The formal inversion of the formulae (34-37) leads to: 

25 *M=(^(m*)-«.£.M-^<--)) , (39) 

= K^—2 L M - (40) 

Z -r €Lj 

W = 4 2 2 (-^L 0 (z) - a.L^z) + L 2 {z)\ quo , (41) 

30 + v - ^ ; 

F:j( = ) = (9+ 4 q2 , 2 (aiL-o(z) + (1 - + a^MJ quo z~ l (42) 

It is then enough to show that the quotients by z A are exact. Since F(z) is monic, L 0 (z = 0) = 
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I, Lj(z = 0) = -a x by defining of a x and according to the relationship (31), L 2 (z - 0) = -a x 12. 
The straight line members of the equalities defining F 2 (z) and F 3 (z) therefore cancel out in z = 
0 and the quotients are exact. It can then be ascertained that the polynomials F r i = 0, 3 
meet the appropriate conditions of symmetry. 
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2 - Equivalent decomposition 

Noting the following matrix as Rf a): 

R(*) = 

and using the formula (32) we get: 

Fo(z) ( 

) PZP T R(c 2 ) PZP T 



cos a — sin a 


0 


0 


sin a cos ck 


0 


0 


0 0 


1 


0 


0 0 


0 


1 



(43) 



FM 



V 1 PJa 



PZP T R(a n ) PZ 



1 

Q-n+l 
1 

0 



15 In the sending diagram, Ma is used to denote the block associated with the right side vector. 



The matrix Z i s then introduced. This matrix is defined by: 

0 0 0 2 

o 



: PZP J 



1 



0 



2 0 taking account of: 



' 1 ' 












1 




0 


0 




0 



(44) 



(45) 



= /%)ZR(flo) Z ... R(a n ) Z 



1 



(46) 



we obtain, with a multiplier factor g, the breakdown: 

. ^(-) . 

Operational complexity 

The equation (46) directly gives the structure of the making of the half-Nyquist sending 
and reception filters. It corresponds to a system with four inputs and four outputs which, in 
the diagrams of Figures 6 and 7, take the place of the set of four polyphase filters. It will be 
noted that, in the case of the sending filter, the outputs have to be inverted with respect to the 
writing (46) which on the contrary corresponds precisely to that of the reception filter. Apart 
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from this difference of detail, the structure for each of these 4n + 2 order filters will therefore 
correspond to: 

• n matrix blocks of the Z and R(a,) types, 

• one P r matrix block, 

• one vector block [1, a n+] *Q, 0] T , 

• a multiplier g which is to be taken into account only for sending or reception. 

Let us now make a more detailed examination of the cost of making each of these 
elements. 

• Matrices RCo^) 

The matrices R(aj), 1 < 1 < n are block matrices which can take the form: 

RM=( o' ?) ( 47 > 

where I is the unity matrix and G Y is the matrix of rotation for the angle a„ namely: 



sin(a/) cos^ 

The circuit arrangement cost is therefore quite simply that of the lattice corresponding to 
the rotation G ; , namely four multiplications and two additions. 

• The matrix Z. 

It can be used for the association, with an input vector e, of the output vector s by the 

relationship s = Ze which is expressed by the following four equations: 

•so = f [(e 0 - e 2 ) + - _1 (e 0 + e 2 )] 

sl= |[(e 2 1 -c 0 ) + -- 1 (e 0 + e 2 )] (49) 



•S3 = e 3 

The structural diagram of a block of this kind is shown in Figure 14. It can be 
specified that the multiplications by 1/2 can be reduced through a simple shift of binary data. In 
addition to these two shifts, the cost is therefore simply four additions and one sign inversion. 
• The matrix P T 
30 It takes the form: 
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/ a 0 

0 1 

a 0 

V 0 0 



0 \ 

0 
0 

1 J 



(50) 



(51) 



with a = -= . This embodiment can also be deduced trom the input/output equations which in 
5 the present case are written as follows: 

s 0 = a(e 0 -f e 2 ) 
•si = e 1 

&2 = a(-e Q + e 2 ) 
S3 ~ e 3 

immediately it is therefore possible to then deduce the diagram of Figure 15. 
1 0 The operational cost is two multiplications, two additions and one sign inversion. 

• The vector block [1, a n+1 , 0, 0] T 

Its circuit arrangement cost is the equivalent of only one single multiplication. 
The full circuit arrangement diagram according to the equation (46) therefore 
corresponds to a system of four inputs and four outputs as shown in Figure 16. Naturally, 

1 5 this system is found at sending and at reception. The order of complexity of each of these 

filters is therefore An MPU and 6n APU. 

Thus, as already stated, at sending it is necessary to invert the outputs. In the drawing 
showing the making of the device in Figure 17, this operation is symbolized by the block J 
corresponding to the antidiagonal matrix, namely: 

20 / 0 0 0 1 \ 

0 0 10 
0 10 0 

V 1 0 0 0 , 

The reception diagram for the 4n + 2 order is shown in Figures 18 and 19. It is 
deduced from the following computations. 

2 5 The matrix 2 i s introduced, defined by : 

|(1 + 0 |(1 0 

0 0 0 z- 1 

iCl-z- 1 ) 0 1(1 + Z- 1 ) 0 

0100 

It is the matrix obtained in taking the mirror elements of the elements of the matrix Z . We 
30 have: 



(52) 



Z = 



z z = 
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where I 4 is the 4th order identity matrix. 

Then the following identities are verified: 

PP T =I 4 , R(a)- 1 =R(a) T = R(-a) 

for any value of a. 

By noting the matrix product as M: 

M = P T R(ci)Z.,.R(a n )Z 
and the matrix product as N : 

N = ZR(-On) • - - ZR(— ai)P 



we therefore have: 



NM = z~ n I 4 , 



and consequently : 



We then introduce the matrix K defined by: 

"10 0 0 



r *i>(*) - 




' 1 _ 


^) 








= gz n 


0 


L F z {z) _ 




0 



K = 



0 10 0 
10 0 0 
.0001. 

Figure 19 gives the drawing for the making of the block corresponding to the matrix N 
2 5 and Figure 20 is the for the making of the reception filter resulting therefrom. If the constant g 
is equal to 1/(1 + a n+1 ), then a signal transmitted by the sending system of Figure 17 and then 
received by the system of Figure 20 produces an identical signal with a delay of n + 2 
samples. 
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APPENDIX 4 



Single-lattice reception diagram 



If we assume a back-to-back operation, it is indeed possible to simplify the diagram of 



5 Figure 10 by using the equation that fulfils a null ISI. The single-lattice diagram of 
Figure 10a is then obtained. In this case, only the null ISI and the phase linearity of the 
sending filter are provided. 

The validity of this drawing of Figure 10a can be verified from the equation (37) of the 
patent after multiplication by r 1 . For a so-called back-to-back system, the input signal of the 
1 0 reception filter is given by: 



Y(z)=z- 1 [F 0 (z*) + z- 1 F 1 (z*) + 



z^F 1 (z 4 ) + 



z- 3 F 0 (z*)]X(z i ) (53) 



Referring to Figure 10a, the signals after decimation are given, from top to bottom, by: 



15 




At input of the direct lattice, the signals are therefore given, from top to bottom, by: 



20 



z-i[F 0 (z)_+ F 1 (z)}X(z) et 



z^[Fo(z)~F 1 (z)]X(z). 



The output signal is therefore: 



S(z) = 2gz- 1 [F Q (z)F Q (z) + F l (z)F 1 (z)}X(z) = 2gz- l [-yz- n ), (55) 



25 



with 



7=n( i +°. 2 )- 



It is thus verified that for g = 1/y, the output signal is identical to the input signal except 
3 0 for the n + 1 sample delay. 
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APPENDIX 5 
The case of the N = 4n + 1 order 

Preliminary observation - To preserve a lattice with n + 1 cells with coefficients 
5 with 0 < i < n, as defined for the 4n + 3 order, we present the case of the 4n + 5 order. 

Let L ; , i = 0, 3 be the polyphase components of a null ISI 4n + 3 order L(z) filter. 
It is possible to associate, with this filter, a lattice built as in Figure 8a of the patent with a set 
of coefficients a t , i = 0, n. We have L 2 (z) - 1 } (z) and L 3 (z) = L 0 (z). 
The filter F(z) whose components F $ i-0,...,3 are defined by 

10 

(56) F 0 (z) = Lo(*), 

(57) F^z) = z- l Lo(z), 
15 (58) W = ^(*)> 

(59) F 3 (z) = iiW, 

is a null ISI, 4n + 5 order linear phase filter. Furthermore, the Theorem 2 in Appendix 2 of 
20 the patent establishes the fact that such filters are necessarily obtained in this way. Of course, 
the synthesis of such a filter is done so as to optimize the coefficients a x , i - 0, n to obtain 
the right frequency specifications for the filter F(z), not for the filter L(z) ! This therefore leads 
to the obtaining of diagrams for sending and reception associated with F(z) from the lattice built 
with the coefficients a t and the associated inverted lattice. 
25 The type I polyphase breakdown applied to this filter is expressed by the following 

expression: 

F{z) = L 0 (z 4 ) + z -*[z-*Lo(z*)] + + ^%(z 4 ) (60) 
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It can be ascertained that, for the sending part, this equation may be expressed by the 
lattice diagram of Figure 20 wherein a delay is applied. This delay is given by z' 3 at output of 
F(z). 



3 1 



For the reception part as in the case of the 4n + 3 order, it is possible to make a 4-lattice 
system that ensures the totality of the properties sought. 

Similarly, simplifications appear if we consider the system known as the back-to-back 
system. The input signal of the reception filter can then be written as follows: 

Y(z) = [z-*L 0 &) + z-^Uz*)] + + z-«Uz*)]X(z*) (61) 

To recover the input signal, apart from a delay, we then implement the device of 
Figure 21. 

Referring to Figure 21, the signals after this decimation are, from top to bottom, given 

by: 

Y(z) 4-4= z-*Lo(z)X(z) t z-'Y(z) U= z-'L 0 (z)X(z), 
z- i Y(z)U=z- 2 L 1 (z)X(z) 1 z^Y(z)U=z-^L 1 (z)X(z) 

From top to bottom, the signals at input of the lattice are respectively given by: 

z-*[Lo(z) + U[z))X(z) et z-*\Lq{z) - L^z)]X{z). 

Thus, as in the case of the 4n + 3 order but this time with a processing delay of n + 2 
samples, we have an output signal which, for g = l/(2y) is identical to that of the input. 
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APPENDIX 6 

1 - The evaluation of the complexity of the reception filters 

Two cases have to be considered depending on whether it is sought or not to keep all 
5 the properties sought. Thus, for the 4n + 3 and 4n + 5 order filters: 

• To keep null ISI, the matched pair, the phase linearity of the sending and reception 
filters, as shown in Figures 10, four lattices are used. These four lattices work at the lowest 
rate. Given the fact that for each lattice two multiplications and additions are necessary, except 
for each first cell (one multiplier and one adder), the operational complexity is (8n + 4) MPU 

1 0 and APU, plus a sign inversion. 

• For a back-to-back system, simplifications are possible but only null ISI and the 
phase linearity of the sender filter are guaranteed. For the simplified 4n + 3 and 4n + 5 order 
reception filters, the operational complexity corresponds then to that of a single lattice plus that 
of the three adders, the multiplier by g and the inverter. The computations performed at the 

1 5 lowest rate can therefore be estimated at (2n + 3) MPU, (2n + 5) APU and one inversion. 

2 - A method of generic synthesis 

As compared with the initially proposed method of synthesis, we now have a synthesis 
method that can be used, for all three types of solutions, to meet a frequency specification such 
20 as the one shown in Figure 2. 

For the filter F(z) with a length L, different from 4n + 1, in other words for an N order 
different from 4n and for a fixed value of the fall-back factor p, we seek to minimize the cost 
function: 

25 = su V {w P (1 - |F(l)j) 2 ,w c {^f - \F{^)\j ,»ssup |F(e^)| 2 }. 

We use the FSQP (feasible sequential quadratic programming) algorithm developed by 
the team led by A.L. Tits [5] to search for the minimum of a set of non-linear even cost 
functions subjected to general, even and non-linear constraints. Since this method is a method 
30 of local optimization, the choice of an initial point is necessary. For a given weight w in the 
stopping band and a weight equal to 1 at the starting point and in the Nyquist frequency, 
namely 7t/8, we calculate an initial filter F imt that is optimal for the norm of the minimal value. 
This filter F mit is then used for direct computation of a set of lattice coefficients for a filter 
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producing a matched pair with null ISI. Since F'"' does not itself produce a matched pair with 
null ISI, the identification between the two sets of coefficients, the transversal coefficients of 
F init and the lattice coefficients is not exactly achieved. We therefore obtain a different filter, 
referenced ¥ imt , whose lattice coefficients are considered to be the initial point of the 
5 optimization problem described by (63). For this problem, we fix w P = 1, w c - 2, w s =0.5. 

Furthermore, for any lattice structure optimized for a given value p 0 , we also use a 

continuation algorithm to obtain a lattice structure optimized for a value p } , that is different in 

taking a sequence of intermediate values for p: the result of the lattice structure optimized for a 

value of p in the sequence is the initial point for the following value of p. 
1 0 Figure 22 shows a higher attenuation in the stopping band than could be obtained for a 

fixed length and fall-back factor. 

Let us consider for example the values L = 43 and p = 0.5: the best attenuation 

obtained by direct optimization is equal to 52.25 dB. The frequency response of the 

corresponding sending filter F°p { is given by Figure 23. The transversal coefficients and the 
1 5 coefficients of the lattice structure are given in Table 3. 
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Transversal coefficients 



5 



fo 


—1.033698 10 4 


In 


8.074275371845 10 3 


fl 


O ^AC A 0701 n A 1 A— 4 

— o.t>yo4z/ zl / 4 iu 


fl2 


1.0/ /4oO/ /4z±y4 1U 


h 


6.60549905248 10~ 4 


fiz 


1.2603068512975 10~ 2 


h 


-1.09979829792 10~ 4 


fu 


-6.358804372375 10~ 3 


h 


4.31563793685 10~ 4 


/l5 


-3.3366501346017 10~ 2 


h 


6.08998879046 10~ 4 


fie 


-4.6538674170018 10~ 2 


fs 


6.39123817678 lO^ 4 


fl7 


-2.1044564354449 10~ 2 


fi 


-7.82626716253 10~ 4 


fis 


5.1214636348857 10~ 2 


fs 


-3.761395765875 10~ 3 


/l9 


1.51186813816588 10" 1 


U 


-4.983816032565 10~ 3 


/20 


2.39358058463701 10" 1 


fio 


-9.24070721312 10~ 4 


M 


2.74720061562686 10" 1 



lil 



Lattice coefficients 



20 



25 





-3.1848708937 


&2 


2.648518499 10" 1 


013 


-2.9385797132 


a 4 


-9.349457966 10" 1 


Q'o 


9.246893837 10" 1 


ore 


-1.308356726 10" 1 


a 7 


3.6683584963 


a s 


-6.126257549 10" 1 


as 


-1.2465410017 


otio 


-1.439677875 10" 1 


a> n 


-3.5749582735 
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Table 3 : Lattice and transversal coefficients of the example with L = 43 and p = 0.5 (f 42< = fj 
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CLAIMS 



1. Method for the making of a digital Nyquist filter with null inter-symbol 
interference designed to process a physical signal transmitted between a sender and 
a receiver through a transmission channel, 

said filter being an Nth order P(z) = F^(z) symmetrical filter implementing an 
oversampling factor M=4 and forming a matched pair comprising a sending filter 
(12) and a reception filter (15) whose polyphase breakdown of F(z) can be written 
as follows: 

F(z) = F 0 (z 4 ) + z~ 1 F 1 (z 4 )+ z' 2 F 2 {z 4 ) + z' 3 F 3 (z 4 ) 

characterized in that N is different from 4n, n being an integer: 
and in that : 



F being the mirror symmetry of F and y being a non-null constant. 

2. Method according to claim 1, characterized in that N is equal to 
4n+3 or 4n+l and: 

said sending filter (12) performs an interpolation (121) by a factor M = 4 

and has a circuit arrangement corresponding to a polyphase breakdown 

known as the type II breakdown, such that: 



If N=4n+2, 



If N=4n+1, 



If N=4n+3, 



F 1 (z)F 1 (z)+ z- 1 F 2 (z)f 2 (z) = TZ _n 
2F 0 (z)F 0 (z) + F 2 (z) +z _1 F;(z) = tz" 
F 0 (z)F 0 (z) + F ] (z)F 1 (z) = Y ^ n 



F F(z) = [z" 3 z -2 z _1 l] 



'F 0 (z 4 )' 
F,(z 4 ) 
F,(z 4 ) 
F 0 (z 4 ) 



and said reception filter (15) performs a decimation (152) by a factor M = 
4 and has a circuit arrangement corresponding to a polyphase breakdown 
known as the type I breakdown, such that: 



F(z)=[F 0 (z 4 ) F,(z 4 ) P^F^z 4 )] 



3. Method according to any of the claims 1 and 2, characterized in that, 
in said sending filter (12), a filtering step followed by a step of interpolation by a 
factor of M=4 is performed. 

4. Method according to any of the claims 1 a 3, characterized in that, in 
said reception filter (15), a step of decimation by a factor M=4 is performed, 
followed by a filtering step. 

5. Method according to any of the claims 1 to 4, characterized in that 
said sending filter (12) and/or said reception filter (15) have a structure in the form 
of at least one lattice. 

6. Method according to claim 5, characterized in that said sending filter 
(12) and said reception filter (15) are each constituted by a pair of polyphase 
components respectively given by the following equations : 



where g is a non-null constant of standardisation and Oj are real coefficients. 

7. Method according to claim 6, characterized in that it implements a 
two4attice structure. 

8. Method according to claim 6, characterized in that it implement a 
single4attice structure working at a double frequency. 

9. Device for the filtering of Nyquist digital signals with null inter- 
symbol interference designed to process a physical signal transmitted between a 
sender and a receiver through a transmission channel, 

based on an Nth order P(z) = F 2 (z) symmetrical filter implementing an 
oversampling factor M = 4 and forming a matched pair comprising a sending filter 
(12) and a reception filter (15), 




o 



A(«) = 



and \z)= 



with : 
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the polyphase breakdown of F(z) of this symmetrical filter being written 
follows: 

F(z) = F 0 (z 4 ) + z-'F^z 4 )* z- 2 F 2 (z 4 )+ Z - 3 F 3 (z 4 ) 

characterized in that N is different from 4n, n being an integer, 

IfN=4n+l, F 1 (z)f 1 (z)+ z- 1 F 2 (z)F 2 (z) = yz" n 
If N=4n+2, 2F 0 (z)F 0 (z) + F~(z) + z- 1 F;(z) = YZ~ n 
IfN=4n+3, F 0 (z)F 0 ( 2 ) + F 1 (z)F 1 (z) = Y z _n 

F being the mirror symmetry of F and y being a non-null constant. 



ABSTRACT 

The invention pertains to a method for the making of a Nyquist digital filter with 
null inter-symbol interference, designed to process a physical signal transmitted 
between a sender and a receiver through a transmission channel, 
said filter being a symmetrical N order P(z) = F^(z) filter implementing an 
oversampling factor M=4, and forming a matched fair comprising a sending filter 
and a reception filter, 

the polyphase breakdown of F(z) being written as : 

F(z) = F 0 (z 4 ) + z" 1 F 1 (z 4 )+ z~ 2 F 2 (z 4 ) + z' 3 F 3 (z 4 ) 
characterized in that N is different from 4n, with n as an integer, 
and in that : 

If N=4n+1, F 1 (z)F 1 (z)+ z- 1 F 2 (z)F 2 (z) = Y z" n 

If N=4n+2, 2F 0 (z)F 0 (z) + F;(z) +Z - 1 F;(z) = Y z" n 

If N=4n + 3, F 0 (z)F 0 (z) + F , (z)F x (z) = fz' n 

F being the mirror symmetry of F and y being a non-null constant. 
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ifjcrcby authorize them to act and rely on instructions from and communicate directly with the pcrson/assignee/attomey/finn/ organization 
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